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Abstract 

A  computational  fluid  dynamic  (CFD)  analysis  was  used  as  a  useful  tool  to  obtain  numerical  relationships  among  similar  flow-field  designs 
having  different  active  areas.  Assuming  that,  at  fixed  current  density,  the  porous  layer  velocity  distribution  calculated  for  the  smallest  size 
geometry  was  the  optimal,  a  scale-up  methodology  has  been  proposed  and  implemented  in  software.  It  allows  acquisition  of  the  scaled-up 
flow-field  by  multiplying  the  reference  cell  parameters  by  suitable  factors  (scale  factors).  A  calibration  methodology  has  also  been  proposed 
by  the  introduction  of  a  feedback  factor  evaluated  by  CFD  analysis.  A  first  application  of  this  procedure  has  been  carried  out  to  scale-up  a 
5-125  cm2  cell. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 

Keywords:  CFD;  Scale-up;  Fuel  cell;  Flow-field  design 


1.  Introduction 

The  polymer  electrolyte  fuel  cell  (PEFC)  is  becoming  an 
interesting  power  source  in  various  application  fields,  from 
portable  units  to  large  stationary  systems.  This  wide  range  in 
power  output  implies  a  great  variation  of  the  active  area  size 
depending  on  the  application.  Generally,  the  scale-up  of  flow- 
field  geometry  from  the  lab  scale  single  cell  to  the  stack  gives 
performance  losses  due  to  the  fluid  dynamic  changes  intro¬ 
duced  by  the  scaling  procedure.  In  the  PEFC,  both  hydrogen 
and  alcohol  fed,  the  interaction  between  the  flow-field  and 
the  electrode  diffusion  layer  plays  a  primary  role  [1-3].  Then 
when  a  cell  configuration  (flow-field  plus  electrode)  has  been 
optimised  in  the  laboratory,  the  maintaining  of  the  reactant 
distribution  and  consequently  of  the  performance  requires  a 
new  optimisation,  increasing  of  the  realisation  time  and  the 
cost  of  the  system. 
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Availing  ourselves  of  computational  fluid  DYNAMIC 
(CFD)  analysis,  a  computer  aided  fuel-cell  design  method¬ 
ology  was  realised  able  to  reduce  the  performance  loss  and 
to  limit  the  number  of  lab- tests.  The  method  involves  the  pa- 
rameterisation  and  the  automatic  generation  algorithm  of  the 
most  common  flow-fields  as  previously  reported  [4]. 

In  the  present  work,  the  algorithm  was  used  as  part  of 
a  computer  aided  selection  method  and  flow-field  scale-up. 
The  design  selection  methodology  is  explained  by  describing 
the  main  program,  the  subroutines  and  their  interaction  sepa¬ 
rately.  The  software  application  to  the  scale-up  from  a  starting 
cell  of  5-125  cm2  is  also  reported.  It  was  observed  that  the 
input  variable  constrains  influence  the  results.  A  method  for 
algorithm  calibration  was  introduced  by  means  of  a  check 
parameter  taking  into  account  the  mass  transfer  towards  the 
electrode.  The  CFD  analysis  has  been  used  to  evaluate  this 
parameter  and  the  methodology  general  aspects. 

2.  Design  selection  method 

An  automatic  scaling-up  design  procedure  demands  the 
definition  of  parameters  describing  the  object  and  their  vari- 
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Nomenclature 

MEA  geometrical  parameters 
AMEA  electrode  surface  area 
LMEA  MEA  width 

Flow-path  parameters 
LCAN  sub-channel  width 
LCOS  rib  width 
HCAN  channel  height 
NDIV  sub-channel  number 

Flow-path  generation  parameters 

Coll  “connection  module”  height 

Hmea  MEA  height 

Lmea  MEA  width  updated  value 

Lser  serpentine  equivalent  width 

Lung  “base  module”  height 

N  reference  value  for  modules  repetitions 


the  cathode  side  of  an  H2/air  cell  was  considered,  and  the  se¬ 
lected  operative  conditions  were:  J=  500  mA  cm-2  (current 
density),  T=80°C  (operating  temperature),  p  =  room  pres¬ 
sure,  RH  =  100%  (air  relative  humidity),  x  =  l  .5/2  (H2/air  sto¬ 
ichiometric  ratio).  According  to  equation  and  symbols  used  in 
reference  [4]  (see  Appendix  A),  the  assumed  operative  condi¬ 
tions  are  used  to  calculate  the  stream  inlet  velocity  as  follows: 


Fmedr 


1  J  AMEA 

- Mr - 

Min  nF 


(2.1) 


where  p  is  the  stream  inlet  density,  A[n  is  the  flow-field  inlet 
area,  Mr  represents  the  reactant  molecular  weight,  F  is  the 
Faraday’s  constant,  n  the  electrons  number  involved  in  the 
semi-reaction  and  AMEA  is  the  total  active  area. 

In  Fig.  2  the  reference  flow-field  is  shown  and  its  charac¬ 
teristics  reported  in  the  table  on  the  right  side,  where  the  nmed 
and  other  fluid  dynamic  parameters  values  are  calculated  at 
the  selected  operative  conditions. 


2.7.  MAIN  routine 


ation  range  (constrains).  In  the  case  of  fuel  cells,  the  fol¬ 
lowing  parameters  related  to  the  channel  gas  flow  path  have 
been  defined:  meander  type  (serpentine,  ribbed,  etc.),  active 
area  surface  and  geometry  (square,  rectangular,  etc.),  section 
(square,  circular,  triangular,  etc.)  width  and  depth  of  the  chan¬ 
nels,  rib  width,  etc.  In  addition  to  geometrical  characteristics, 
fluid  dynamic  properties  have  to  be  also  considered,  because 
a  correct  reactant  distribution  over  the  electrode  active  area 
is  fundamental  to  reach  a  good  fuel  cell  performance  [1-3]. 
Starting  from  the  above  mentioned  parameters  the  computer 
aided  scale  up  method  permits  us  to  obtain  in  the  active  area 
a  fixed  range  of  the  flow-field  giving  performance  as  close  as 
possible  to  the  reference. 

The  developed  software  consists  of  four  routines  (Fig.  1), 
the  main,  the  flow-field  generation  software  (FGS),  the  selec¬ 
tion  algorithm  software  (SAS)  and  the  post-graphic  output 
(PGO). 

The  MAIN  routine  elaborates  the  active  area  within  the 
selected  range  and  operative  conditions  inputs  and  sends  the 
calculated  construction  parameters  to  the  FGS  that,  accord¬ 
ing  to  the  methodology  explained  in  [4],  generates  all  the 
possible  serpentine  like  flow-fields:  single  (SS),  ribbed  (RS) 
and  multi- serpentine  (MS).  Then,  the  SAS  subroutine  selects 
the  design  which  matches  geometrical  and  fluid  dynamics 
constraints.  At  each  ith  step  the  selected  flow-field  is  stored 
in  a  database  and  successively  used  for  diagram  generation 
by  the  post-processing  graphical  output  (PGO)  routine.  PGO 
also  generates  the  scale  factor  and  diagrams  as  a  function  of 
the  cell  active  area. 

A  conventional  single  serpentine  flow-pattern  was  se¬ 
lected  as  the  sample  case  because  it  is  the  most  common  flow- 
field  in  the  PEFC.  In  particular,  the  commercially  available 
5  cm2  active  area  used  for  the  lab  test  cell  was  a  single  serpen¬ 
tine  with  a  rectangular  channel  of  1  mm  in  width  and  1  mm  in 
depth  and  10  180°  turns  (11  passages  total).  As  a  sample  case 


The  MAIN  routine  spans  the  active  area  from  AME Amin  to 
AMEAmax  with  a  defined  increment,  A  AMEA.  MAIN  input 
variables  are  the  cell  operative  conditions,  the  total  MEA 
active  area  range  and  the  step  to  be  used  (Table  1).  For  each 
AMEA t  the  routine  calculates  all  the  variables  to  be  send  to 
the  FGS  routine. 


AMEAmi„<  AMEAj  <  AMEA 
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Fig.  1.  Scale-up  software  flowchart. 
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Input  parameters 


MEA  geometrical  parameters 

AMEA 

5  cm2 

LMEA 

22  cm 

Flow-Path  parameters 

LCAN 

1.0  mm 

HCAN 

1.0  mm 

LCOS 

1.0  mm 

NDIV 

1 

Output  parameters 


Flow-Path  generation 
parameters 

Design  coefficients 

% 

271.90  mm 

Coll 

1.0mm 

c 

0.88 

Lser 

1.0mm 

X 

0.54 

Lung 

21.81  mm 

Fluid  dynamic  quantities 

N 

11 

Lmea 

21.0  mm 

V me <J 

3.32  m/s 

Dh 

1 .00  mm 

Hmea 

23.81  mm 

Re 

82.54 

Ap 

9.60  mbar 

Fig.  2.  Sketch  of  5  cm2  reference  cell  and  tables  of  its  parameter  values. 


The  inputs  are  used  to  calculate  the  reactant  stream  mass 
flow  rate  and  its  thermodynamic  properties,  namely  fluid  den¬ 
sity  and  fluid  viscosity  and  ith  reference  inlet  area: 


A 


inrcf(/) 


rhi 

PVme  dref 


(2.2) 


the  generation  of  all  the  possible  designs  having  geometric 
characteristics  inside  the  defined  ranges.  The  FGS  takes  into 
account  the  outputs  coming  from  the  MAIN  (see  Table  2) 
and  the  geometric  constrains  to  calculate  the  solution  matrix 
dimensions  by  the  following  equations: 


where  pmedref  represents  the  average  inlet  velocity  of  the  ref¬ 
erence  cell,  therefore  Ainref(/)  stands  for  the  channels  cross- 
sectional  area  that  gives  an  inlet  velocity  value  equal  to  that 
of  the  reference  cell. 

The  number  of  iterations  depends  on  AMEAmin, 
AMEAmax  and  AAMEA  as  follows: 


/  AMEAmax  AMEAmin  \  , 

W cycles  -  AAMEA  ) 


(2.3) 


Fluid  density  and  viscosity  were  calculated  by  simple 
mathematic  media  assuming  that  the  stream  is  a  mixture  of 
two  perfect  gas,  air  and  water  vapour. 


2.2.  Flow -field  generation  software  FGS 

The  methodology  for  automatic  serpentine  like  flow-field 
generation  previously  evolved  [4],  was  implemented  to  allow 

Table  1 


Main  program  input  and  output  variables 


Input 

Output 

Total  active  area  max  value,  AMEAmax 

ith  total  active  area,  AMEA, 

Total  active  area  min  value,  AMEAmin 
Total  active  area  step  value,  AAMEA 

Ref.  inlet  area,  Ajnref(0 

Operative  conditions 

Current  density,  J 

Gas  flow  rate,  m; 

Operating  temperature,  T 

Fluid  density,  p 

Back-pressure,  p 

Relative  humidity,  RH 

Stoichiometric  ratio,  x 

Fluid  viscosity,  n 

nr  =  int 


A 


in¬ 


max 


AAin. 


(2.4a) 


NDIV 


nc  = 


max 


ANDIV 


(2.4b) 


Table  2 

FGS  parameters 

Input  Output 


Active  area  parameters 

Electrode  surface  area,  AMEA ; 
Ref.  ith  inlet  area,  Ainref(/) 

Inlet  steam  characteristics 

Gas  flow  rate,  m 
Fluid  density,  p 
Fluid  viscosity,  /i 

Flow-path  parameters 

Channel  height  (fixed),  HCAN 
Rib  width  (fixed),  FCOS 


Flow-path  parameters 

Sub-channel  width,  LCAN, 
Sub-channel  number,  NDIV, 
Flow-path  generation  parameters 
(FGPs) 

ME  A  width,  Fmea  , 

MEA  height,  Hmea, 

Reference  value  for  modules 
repetitions,  A) 

“Connection  module”  height, 

Coif 

Serpentine  equivalent  width,  Lser, 
“Base  module”  height,  Fung, 
Design  coefficients  (DCs) 

“Unrolled”  serpentine  length,  ft 
MEA  aspect  ratio,  0 
Covering  factor,  xi 
Fluid  dynamic  quantities  (FDQ) 
Average  cannel  velocity,  umed 
Hydraulic  diameter,  D h 
Reynolds  number,  Rei 
Total  serpentine  pressure  drop, 

A  pi 
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where  Ainmax  and  NDIVmax  are  evaluated  as  follows: 


^  inmax  _ 

l<A]nrcf(i)  with  0  <  k  <  1 

(2.5) 

NDIVmax  —  int(Ajnmax)  3“  ^ 

(2.6) 

where 

5  =  « 

r 

0 

1 

ifdec(Ainmax)  <  0.5 
ifdec(Aillmax)  >  0.5 

Eqs.  (2.4a)  and  (2.4b)  represent  the  rows  and  columns  number 
of  the  solution  matrix,  respectively,  whereas  Eq.  (2.5)  gives 
the  maximum  value  of  cross-sectional  inlet  area  as  a  multiple 
of  Ainref(/),  and  Eq.  (2.6)  the  number  of  channel  subdivisions 
(NDIVmax). 

Consequently  the  number  of  possible  designs  is: 

^sol  =  nYnc  (2.7) 


be  set  to  avoid  important  pressure  losses  in  the  flow-field. 
Therefore,  ith  sub-channel  width  is  admitted  only  if: 

LCANmin  <  LCAN/  <  LCANmax  (2.1 1) 

A  sub-channel  width  proportional  to  0. 1  mm  was  fixed  for 
manufacturing  reasons,  so  that  the  SAS  algorithm  eliminates 
all  the  solutions  which  do  not  respect  the  subsequent  condi¬ 
tion: 

LCAN/ocO.l  (2.12) 

As  is  well  known,  the  electrode  contact  area  must  be  max¬ 
imised  to  minimise  the  contact  ohmic  losses,  otherwise  a  free 
area  needs  to  assure  reactants  of  diffusion  through  the  elec¬ 
trode  backing  toward  the  active  catalyst  sites.  This  parameter 
is  taken  into  account  by  the  covering  factor  /  which  is  defined 
as  the  free  area  to  total  active  area  ratio,  thus  for  ith  solution 
Xi  value  is  limited  by: 


To  reduce  the  solution  matrix  dimension  the  channel  height 
(HCAN)  and  rib  width  (LCOS)  were  considered  as  fixed 
values,  while  the  sub-channels  number  (NDIV/)  and  width 
(LCAN/)  vary  to  cover  the  ith  inlet  area: 

NDIV/  =  NDIVmin  +  j  A  NDIV,  j  =  1  -nc, 

Am i  —  Ainmin  T  /AAin?.,  /  =  1  —nv  (2.8) 

where  NDIVmin  is  the  minimum  allowed  value  fixed  (=1  for 
simple  serpentine),  A  NDIV  represents  the  increment  (in  our 
case  =  1),  Ainmin  is  a  percentage  of  Ainref(/)  and  AAin.  repre¬ 
sents  the  channel  cross-section  inlet  area  increment. 

Given  Ain/  and  NDIV/,  LCAN /  and  the  other  output  pa¬ 
rameters  are  immediately  determined  for  each  iterative  cycle 
as: 


LCAN  / 


^in, 

NDIV,  HCAN 


(2.9) 


Thus  the  software  for  each  Ain.  that  is  Ain  .  <  A{n.  <  A^. 

xxxi  AAAmin  —  xxxi  —  xxxuidx 

calculates  the  flow  path  generation  parameters  (FGPs)  and  the 
design  coefficient  (DCs)  by  varying  NDIV/. 

The  inputs  and  outputs  of  FGS  are  listed  in  Table  2.  FGPs, 
DCs  and  FDQ  were  calculated  by  using  equations  described 
in  [4]  and  listed  in  Appendix  A. 


Xmin  A  Xi  —  Xmax  (2.13) 

Studies  of  flow-field  and  electrode  diffusion  layer  interac¬ 
tions  [1,2]  suggest  that  the  channels  gas  shorting  depends 
on  the  total  flow-field  pressure  drop  which  is  related  to  the 
serpentine  turns  number,  N.  In  fact,  the  N  parameter  influ¬ 
ences  the  pressure  driven  mass  flow  rate  through  the  porous 
layer,  therefore,  the  SAS  imposes  that  the  N  parameter  for  ith 
generated  flow-field  is  in  the  range: 

N min  A  A/  A  N max  (2.14) 

where  Vmin  and  NmSLX  values  are  empirically  fixed. 

A  further  constraint  is  imposed  from  the  MEA  aspect  ratio, 
f,  ranging  between  a  £min  and  a  £max  to  have  the  desired 
active  area  shape.  For  the  sample  case  the  aspect  ratio  of  the 
reference  cell  was  chosen. 

The  SGS  provides  also  a  selection  of  the  flow-field  type 
and,  for  the  reported  sample  case  the  following  selection 
equations  were  used: 

NDIV  =  1  (simple  serpentine)  (2.15) 

NDIV  >  1  and  Coll  >  Collmax 

(multiple  serpentine)  (2.16) 


2. 3.  Selection  software  ( SAS ) 

Even  if  the  solutions  coming  from  the  FGS  are  mathe¬ 
matically  correct  not  everything  can  be  practically  realised, 
so  that  the  SAS  provides  a  solution  filtering  on  the  basis  of 
geometrical  constraints. 

The  first  constraint  considers  the  maximum  usable  sub¬ 
channel  width,  in  fact  the  MEA  flexibility  causes  a  penetra¬ 
tion  of  the  electrode  into  the  channels  obstructing  the  gas 
flow  path.  So  that,  if  the  mechanical  properties  of  the  MEA 
are  known  [5],  the  suitable  sub-channel  width  could  be  fixed, 
moreover  a  minimum  value  for  the  sub-channel  must  also 


NDIV  >  1  and  Coll  <  Collmax 

(ribbed  serpentine)  (2.17) 

in  fact  the  NDIV  parameter  value  is  a  switch  between  SS 
and  MS/RS,  while  the  Coll  parameter  value  allows  to  choose 
between  MS  and  RS  flow-field,  the  Collmax  can  be  fixed  by 
the  user,  considering  the  mechanical  MEA  properties  and 
flow-path  fluid-dynamic  aspects. 

At  this  point  the  fluid  dynamic  conditions  can  be  imposed 
to  select  a  unique  flow-field  design.  In  the  sample  case  a 
pressure  drop  as  close  as  possible  to  the  reference  flow-field 
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pressure  drop  was  required.  A  pressure  drop  difference  per¬ 
centage  (A p%)  term  representing  the  variation  in  pressure 
drop  from  the  reference  value  was  introduced: 


A  pi  A/?rcf  \ 


A  pi 


■ 


x  100 


(2.18) 


where  the  ith  component  of  the  A p%  array  is  estimated  by 
Eq.  (2.18). 

The  final  solution  is  the  closest  to  the  reference: 


aPsol  =  min(zl  p%) 


(2.19) 


2.4.  The  PGO  software 


0  25  50  75  100  125 


AMEAj  [cm2] 


The  PGO  software  is  used  to  post-process  the  obtained 
solution  matrix,  each  solution  was  characterized  by  seven 
dimensionless  scale  factors. 

Among  the  FGS  parameters  listed  in  Table  1  those  con¬ 
sidered  most  significant  were  chosen  and  the  ratio  with  the 
corresponding  reference  value  was  defined: 


LCAN, 

a  i  = 

LCANref 

(2.20) 

NDIV/ 

Oil  —  - 

-  NDIVref 

(2.21) 

o  Si 

ft  = 

(2.22) 

fo  =  — 

Xref 

(2.23) 

A  pi 

n  =  A 

A/?ref 

(2.24) 

Fmed/ 

Y2  = 

Fmedref 

(2.25) 

Ret 

Yi  = 

ReKi 

(2.26) 

The  a/,  Pi,  Yi,  factors  are  related  to  the  serpentine  geometry 
(in  Eqs.  (2.20)  and  (2.21)),  to  the  MEA  DCs  (in  Eqs.  (2.22) 
and  (2.23))  and  to  the  FDQ  (in  Eqs.  (2.24)-(2.26))  respec¬ 
tively.  The  scale  factors  value  is  the  unit  when  ith  solution  and 
its  reference  correspondent  parameters  are  equal,  whereas  the 
minimum  and  maximum  values  are  strongly  dependent  on  the 
assigned  constrains. 

The  graphic  output  of  PGO  gives  an  immediate  outlook  of 
scale  factor  variation  in  respect  of  active  surface  area.  After 
constrain  calibration,  the  scale  factor  graphics  furnish  all  the 
flow-field  scale-up  parameters  in  the  suited  active  area  range. 

2.5.  Algorithm  application 

It  is  noticed  that  the  proposed  methodology  needs  to  be  cal¬ 
ibrated  because  of  the  output  parameters  (the  scale  factors), 
coming  from  the  PGO,  are  strongly  related  to  the  assigned 
constrain  values.  Thus,  to  select  a  suitable  geometrical  and 


Fig.  3.  f  behaviour  vs.  cell  active  area  for  different  N  range  values. 

fluid  dynamic  constrains  range,  a  feedback  parameter  has  to 
be  introduced.  In  this  work,  to  apply  the  algorithm  to  a  sample 
case,  a  parameter  xjr  was  introduced  and  evaluated  by  CFD 
analysis  according  to  reference  [2] : 

rae  i 

x/f  =  —  x  100  (2.27) 

rhtot 

Eq.  (2.27)  represents  the  percentage  ratio  between  the  reac¬ 
tant  mass  flow  rate  through  the  electrode  (me i)  and  the  total 
mass  flow  rate  (ratot)- 

This  parameter  was  assumed  as  indicative  of  fuel  cell  per¬ 
formance  and  for  the  reference  cell  we  found  xfr  =  0.58. 

For  each  selected  flow-field  corresponding  to  the  ith  ac¬ 
tive  area,  the  i/q  value  was  calculated  and  plotted  versus  active 
area  to  monitor  the  selection  procedure  response  varying  the 
constrains  range  as  shown  in  Fig.  3.  Three  different  runs  were 
performed  changing  the  allowed  module  repetition  number, 
N,  in  the  following  ranges:  7-11,9-13  and  11-16.  The  flow- 
field  with  a  x/f  close  to  that  one  of  the  reference  cell  was 
considered  the  best  solution.  As  it  can  be  seen  the  method, 
with  the  imposed  constrains,  is  highly  sensitive  to  parameter 

N.  Increasing  the  minimum  N  flow-fields  with  greater  pres¬ 
sure  drop  are  selected  and  the  xf  value  increases  due  to  the 
transition  from  RS  to  MS.  Then,  to  have  a  x[s  value  near  to 

O. 58  by  changing  only  the  A  constrains  we  must  select  a  low 
N  for  great  area  and  high  N  for  small  area  cell. 

3.  Results  and  discussion 

The  method  was  applied  to  scale  up  a  serpentine  flow-field 
from  5  cm2  (lab  reference  cell,  Fig.  2)  to  a  125  cm2  active  area 
cell  using  different  flow-field  type  (SS,  RS  and  MS,  Fig.  4). 

The  used  MAIN  inputs  are  listed  in  Table  3,  while  the 
channel  depth  (HCAN)  and  the  rib  width  (FCOS)  were  fixed 
to  1  mm. 

For  the  ith  AMEA  size  the  maximum  inlet  area  A\„v, 
was  set  according  to  Eq.  (2.5)  where  k  was  fixed  equal  to  1, 
the  Ajnref(/)  value  is  calculated  by  MAIN  using  Eq.  (2.2)  and 
sent  to  the  FGS. 
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(a)  (b)  (c) 

Fig.  4.  Sketch  of  serpentine  flow-field  types:  (a)  simple  serpentine,  (b)  n -ribbed- serpentine,  (c)  n-multiple  serpentine,  f  behaviour  vs.  cell  active  area  for 
different  N  range  values. 


The  SGS  selection  criterion  was  assigned  as  follows  (see 
Eqs.  (2.1 1)— (2. 17): 

1.0  <  LCAN*  <  2.5,  0.50  <  Xi  <  0.65, 

N  =  11,  13,  16,  0.88  <  &  <  1.19, 

Coll  <  3(RS),  Coll  >  3 (MS) 

Fig.  5(a)-(c)  shows  the  obtained  a,  /3,  y  scale  factor  ver¬ 
sus  active  area  diagrams,  they  give  an  immediate  idea  about 
the  changing  of  the  FGPs,  DCs  and  FDQ  respect  to  the  refer¬ 
ence  factors  value,  namely  the  unit.  In  fact,  fixing  the  desired 
scaled-up  AMEA  size  the  diagrams  furnish  all  the  informa¬ 
tion  needed  for  flow-field  construction  and  the  main  fluid- 
dynamic  and  active  area  design  characteristics. 

Fig.  5(a)  shows  that  simple  serpentine  has  to  be  taken 
into  consideration  for  small  size  active  area  cells  only  (up  to 
15  cm2).  In  fact,  increasing  the  active  area  size,  the  flow  rate 
proportionally  raises,  but  the  limitation  of  maximum  channel 
width  should  cause  a  dramatic  pressure  drop  increase.  So  that, 
a  serpentine  channel  sub-division  is  necessary.  The  possible 
solutions  for  higher  active  area  are  the  RS  flow-field  from  15 
up  to  26  cm2,  and  the  MS  flow-field  up  to  125  cm2.  Moreover, 
Fig.  5(b)  and  (c)  shows  the  MEA  DCs  and  FDQ  behaviour, 
respectively. 

In  Table  4,  the  scale  factors  for  a  100  cm2  active  area 
cell,  derived  from  Fig.  5(a)-(c),  are  reported.  Multiplying 
these  scale  factors  by  the  correspondent  reference  parame¬ 
ters  (Table  5)  the  scaled-up  cell  parameter  values  (100  cm2) 
were  obtained  (Table  6).  The  residual  parameters  that  are 


Table  3 

MAIN  input  parameters 


Input  parameters 

AMEAmax  (mm2)  125 

AMEAmin  (mm2)  6 

A  AMEA  (mm2)  1 

Operative  conditions  parameters 
/(mAcm-2)  500 

T  (°C)  80 

p  (bar)  Room 

RH  (%)  100 

1.5/2 


Fig.  5.  Scale  factor  profiles  as  a  function  of  the  active  area  cell:  (a)  a\,  a2‘, 
(b)  fi,  p2;  (c)  y i,  y2. 
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Table  4 

100  cm2  active  area  cell  scale  factors 


a 

1.50 

012 

4.00 

Pi 

1.35 

Pi 

1.125 

Yi 

8.677 

V2 

3.333 

V3 

4.000 

Table  5 


reference  5  cm2  parameter  values 


LCAN  (mm) 

1 

NDIV 

1 

S 

0.88 

X 

0.54 

A p  (mbar) 

9.59 

Vmed  (ms-1) 

3.32 

Re 

82.54 

Table  6 


100  cm2  parameter  values 


LCAN  (mm) 

1.50 

NDIV 

4 

S 

1.19 

X 

0.61 

A p  (mbar) 

83.21 

t’med  (rn  S  ) 

11.06 

Re 

330.16 

Appendix  A 

•  Flow-path  generation  parameters 

Coll  =  NDIVLCAN  +  5(NDIV  -  l)LCOS, 
5  =  1  for  MS;  5  =  OforRS 


Lser  =  NDIVLCAN  +  (NDIV  -  l)LCOS 


Lung  =  Hmea  —  2Coll 

/  LMEA  +  LCOS  \ 

N  =  int  -  +  5 

v  Lser  +  LCOS  ) 


Lmea  =  (V  NDIV)  LCAN  +  (VNDIV  -  l)LCOS 


Hmea  = 


AMEA 


Lmea 
Design  coefficients 

A  lib 


§  = 


NDIV  LC  AN 


f  = 


Lmea 

Hmea 


X  = 


A  lib 


AMEA 


Abb 


Acent  T  A coii  +  Ain?out 


not  listed  in  the  tables  can  be  easily  calculated  by  equations  Acent  —  N  NDIV  LCAN  Lung 

reported  in  Appendix  A. 

Ain, out  =  2(LC AN  NDIV  Coll) 


4.  Conclusions 


A  coll  —  2 


(NDIV  LCAN)2  +  (NDIV 


1)  NDIV 


An  automated  computer  aided  method  to  scale  up  a  fuel 
cell  flow-field  was  developed.  This  method  utilised  an  al¬ 
gorithm  to  generate  possible  flow-field  designs  on  the  basis 
of  the  imposed  geometrical  constrains  and  to  select  that  one 
matching  the  fluid  dynamic  target.  To  verify  the  validity  of 
the  method  it  was  applied  to  a  sample  case  by  scaling  up  a 
5  cm2  active  area  lab  cell  up  to  125  cm2,  and  the  generated 
flow-field  designs  were  analysed  by  computational  fluid  dy¬ 
namics.  For  the  selection  of  a  single  solution  the  parameter 
ijr  [2]  that  accounts  for  flow  distribution  was  used  as  a  target 
function,  although  other  selection  criteria  could  be  consid¬ 
ered.  Analysis  of  obtained  results  versus  constraint  variation 
indicates  that  the  evolved  computer  aided  flow-field  design 
method  is  a  useful  tool.  Further  work  is  necessary  to  calibrate 
the  criteria.  Work  is  in  progress  for  experimental  validation 
of  the  obtained  results  and  for  the  extension  of  the  method  to 
other  flow-field  types  and  geometry. 


x  (LCAN  LCOS)  + 
x  (N  -  1) 

Fluid  dynamic  quantities 


/NDIV  LCAN  LCOS 


V 


m 


m 


Lriccl  — 


p Ain  yo(LC AN  HC AN  NDIV) 


4  Across  LCAN  HCAN 

Dh  =  - =  2 


Fch 


LCAN  +  HCAN 


Re  = 


pVmed  D\\ 
fJ 


2m 


\i  NDIV (LCAN  +  HCAN) 


9  a  m 

Apsc  =  X  Re2  ^  =  X 


LCAN  +  HCAN 
2 pD2  *  _  "4p  (LCAN  HCAN)3  NDIV2 
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Apcm  = 

Aptot  =  ^Psc  H-  ^Pc m 
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